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[ We investigate the ground-state electronic properties of the symmetricahy- 

, doped semiconductor quantum well in the presence of a homogeneous electric 

00 ' 

I field. In this paper we examined the effect of the electric field and carrier 

' interaction on the subband structure as a function of the field strength and 

^ , carrier concentration. The many-body effects are evaluated using a local 

' density functional exchange-correlation potential. We find that the electron 

o 
o 



subband energy is reduced as the magnitude of the electric field is increased, 
but it is increased as the surface carrier density is increased. However, the 



H ' separation of the electron subband energies is reduced for the increase in both 

the electric field and surface carrier density. On the other hand, the energy 
separation of the hole subbands is increased as the carrier density and field 
strength are increased. Effect of the exchange-correlation potential on the 
subband structure is found negligibly small in this calculation. The subband 
energies are reduced slightly, increasing their separations little in the presence 
of the local exchange-correlation potential. 
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1. INTRODUCTION 



In recent years extensive efforts have been devoted to investigating the electronic proper- 
ties of semiconductor quantum well structures These studies were motivated by the 
improvements in the crystal growth techniques, such as molecular beam epitaxy and metal 
organic chemical vapor deposition. These growth methods are capable of producing ul- 
trathin layers of semiconductor heterostructures with sharp interfaces of high quality. A 
high-mobility semiconductor quantum well structures can be realized with a use of the 
modulation-doping technique 0. In modulation-doping technique, for example, only the 
AlGaAs layers are doped with impurities and the GaAs layers are undoped. Doped impuri- 
ties (donor impurities in our case) in the AlGaAs layers are ionized and electrons transfer to 
the GaAs layers. In this case one needs to examine the subband structures self-consistently 
taking into account the effects of charge transfer from AlGaAs to GaAs layers and of the 
band bending. 

In the presence of an external d. c. electric field, new structures and phenomena were 
predicted and observed experimentally . In particular, the tunability of the optical prop- 
erties of the structure in the presence of an external electric field were studied considerably 
in detail Far infrared intraband absorption was also observed in GaAs/AlGaAs 

layers with high detectivity approaching that of HgCdTe The quantum well high- 



electron-mobility transistor has also been studied Most of these studies involve some 
form of modulation of the quantum well eigenenergies and eigenf unctions in the presence of 
an electric field for their principal mode of operation. 

Previous investigations have been focused in understanding the effects of the electric field 
and carrier interaction separately employing methods of either variational or perturbation 
calculations in a frame of the single-particle picture. The variational method, although 
computationally the simplest, has drawback of not knowing the accuracy and region of 
validity [jlT[. Bloss investigated only the electric field effects by solving differential equation 
numerically and neglecting the effects of carrier interaction [|12[ . Ando investigated only the 
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effect of carrier interaction in the absence of the electric field [0. However, in practical 
devices, one should take both of these two effects, self-consistently, into account on an equal 
footing. 

In this paper, we investigate the effects of the carrier interaction and electric field on 
the electronic properties of symmetrically-doped GaAs/AlGaAs single quantum well (SQW) 
structures. We evaluate the electron and hole subband structures self-consistently within 
a Hartree approximation in the presence of a d. c. external electric field perpendicular to 
the interfaces. Also, a simplified local density functional exchange-correlation potential is 
employed to include the many-body effects such as exchange and correlation interactions in 
the calculation. In section 2, we describe a self-consistent numerical calculation of subband 
structures of the SQW in the presence of a d. c. electric field and carrier intercation. In 
section 3, we present the numerical results of GaAs/AlGaAs system. The subband energies 
are presented as a function of surface carrier density for several different values of the field 
strength. Finally, we conclude the work in section 4. 



2. THEORY 



We consider a symmetrically-doped GaAs/AlGaAs SQW structure in the presence of an 



external d. c. electric field within the effective mass approximation |T^. Here we assume that 
the effective mass encounters the effect of bulk band structure and we ignore the mismatch 
of effective masses in the well and barrier regions. The motion of carriers (electrons or holes) 
in the z direction, which is the growth direction of the SQW and is perpendicular to the 
interface, is quantized due to the one-dimensional confining potential well. And, hence, the 
carriers possess only two degrees of freedom along the x and y directions. Therefore, the 
single-particle wave function describing the motion in the xy-plane is simply plane wave like 
and the wave function in the z direction is the solution of a one-dimensional Schrodinger 
equation 

. d^ 



2m* dz'^ 



^Vt{z)\i,{z)=Ei,{z). (1) 



Here m* is the effective mass of a carrier in the quantum well and the effective potential 
energy function Vt{z) subjected to a single particle in the system is given by 

Vt{z) = V^{z) + Vf{z) + Vs.c.{z). (2) 

In equation (|^) V^^(2;) is the rectangular well potential energy associated with the band-gap 
mismatch at heterointerfaces and Vf{z) is the linear potential energy due to the constant 
electric field applied along the growth direction. And Vs.d^) is the energy of a given carrier 
due to the interactions with ionized impurities doped in the barrier material and with the 
other remaining carriers in the quantum well. The interaction potential energy Vs.d^) 
should be determined self-consistently. 

If one simply replaces the wide band-gap material (AlGaAs layer) with a simple potential 
barrier of height Vo, Vw{z) is written by [see figure 1]: 

I 0, \z\<L/2 
K(^) = (3) 
[ Vo, \z\ > L/2 

where L is the width of the gap region (GaAs layer) i.e., the width of the quantum well. 
In this work the band-gap mismatch Vo is assumed to be 60% (or 40%) of the band-gap 
mismatch of AlGaAs and GaAs for the conduction (or valence) band [|l^. In the presence 
of an external d.c. electric field F in the z direction, the motion of the electrons are further 
confined in addition from that in the absence of the field. The potential energy due to the 
electric field is given by 

Vf{z) = TeFz (4) 

where — is for an electron and + is for a hole. In this paper the electronic charge is — e 
(< 0). As one increases F, carriers confined initially by Vw{z) would tunnel out of the 
well and therefore lowering their potential energy. If, however, the electric field is not 
extremely strong, the quantum-confined states would have a long lifetime and can therefore 



be considered as quasi-bound |Tl . 
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In modulation-doped SQW structures of GaAs/AlGaAs, electrons are transferred from 
the impurities doped AlGaAs to the GaAs region so as to lower the potential energy of 
electrons. As electrons transfer to the well region, the conduction and valence band edges 
are modified from the fiat ones. This band-bending effect is to be included in the calculation 
by solving the Schrodinger equation with a use of a self-consistent potential energy Vs,c.{^)j 

VsA^) = Vh{z) + VUz)- (5) 

Here Vh{z) is the Hartree interaction energy, which is the sum of two contributions; the 
electrostatic interaction of a given carrier with the other carriers in the quantum well and 
that with ionized impurities remaining in the barrier region. Vxc{z) is a single particle 
version of the exchange-correlation energy, which is not included in Vh{z). In the Hartree 
approximation one only considers Vh{z) and neglects V^dz) in Vg.d^)- The Hartree term 
Vh{z) is a solution of a Poisson equation. In this work, we consider the case that impurities 
are doped symmetrically with a bulk density A^^^ in each barrier region on both sides of a 
SQW [see figure 1]. Therefore, the Poisson equation is written by ^ 



d^Vniz) ^[n{z) - N^], L/2 < \z\ < L/2 + d. 



d;s2 



t{ (6) 



^^n{z)^ otherwise 
where the upper (lower) sign refers to electrons (holes). In equation (^ n{z) is the carrier 
density of the SQW and rff, is the width of the impurity-doped region in the barrier. The 
total effective width of the SQW is = L + 2db [see figure 1]. The static dielectric constant of 
the quantum-well material e is assumed to be the same in both the barrier and well region. 
The average two-dimensional {2D) carrier density Ng of the well is given by Ng = 2Nddb. In 
this paper, we confine ourselves to the ground-state, zero-temperature limit of the problem. 
Hence the volume density distribution n{z) is written by 

n{z) = J2Ns,\U^)\' (7) 

i 

where Nsi is the 2D density of carriers in the i^^ subband (i.e., J2i = Ng). If the electrons 
occupy only the ground subband {i = 1), Ng = Ngi [which is the case of weakly-doped 
samples in experiment], and we can rewrite equation (|^) by 



n{z)=N,\M^)\^- (8) 

The Hartree approximation is, in general, known to overestimate the Coulomb repulsive 
force of other electrons in the short range, and, in real system, many-body effects such 
as the exchange-correlation interaction could become important in the study of subband 
structures or intersubband resonance absorption ||T5[. We include the exchange-correlation 
effects in our calculation by introducing a simplified local exchange-correlation potential 
yxc{z) = Vxc[n{z)\ = fixc['n{z)] within the local density-functional approximation (LDA) 
|T6HT8|. Here ^xc is the chemical potential of an inhomogeneous electron gas, which is 



obtained from that of a homogeneous electron gas by replacing the uniform electron density 
by the local electron density n{z) of the inhomogenous system. Several different forms of 
Vxc{z) are proposed by a number of people [0. In our work, we use a simple form suggested 



by Gunnarson and Lundqvist 



2 114 m*e'^ 

VxMz)] = [1 + 0.0545r,ln(l + (9) 

Tiars Ts 2e^h 

where a = (4/97r)^/'^ and is defined in term of n{z) by n{z) = [^(ab'^s)^]^^ with a*^ = 
eh^ /m*e^^ the effective Bohr radius. Therefore, the Schrodinger equation for a single particle 
of the potential energy K.c.(-z) given by equations and becomes the 'so-called' Kohn- 
Sham equation. 

In actual numerical calculation, one first needs to introduce impenetrable potential bar- 
riers at the positions sufficiently far from both sides of the SQW in such a way that the 
positions of the barriers do not affect the solutions of the problem. In our calculation, we 
assume these barriers of infinite height at x = iL^ (~ ±2.5L) ||12[ where —L/2 and L/2 
are the edges of the finite SQW [see figure 1]. We, then, write, within the centeral differ- 
ence approximation ||21| , the second-order differential equation, equation (|I]) by a difference 



equation of unknown il){zi) for each point Zj, a point of the i^^ subdivision of the interval 
[—Lbi Lb]- Therefore, we replace the ordinary differential equation, equation (|I]) by a finite 
difference equation; 
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Z^[ + Vt{zi)4j{Zi)] = E}_^^{zi). (10) 

i=i ^ i=i 

In equation A is the grid spacing, which is the length of each subinterval and n is the 
number of subintervals. The boundary conditions require ipo{zo) = ipn+iizn+i) = 0, where 
Zn = Zo + nA. Writing equation (|T0|) in a matrix form, we get an n by tridiagonal matrix 
eigenvalue equation. This eigenvalue equation is solved by diagonalizing the matrix. In 
this work, we solve this problem self-consistently with an accuracy such that the difference 
in the Hartree potential energies of the i^^ and (z + i)*^ iterations is less than 10~^meV. 
The matrix-size dependence of the eigenvalues is given in figure 2, as an example. It shows 
the electron subband energies in the absence of the field within a single-particle scheme, 
that is, Ns = 0; the solution of equation (|T]) with Vt{z) = Vw{z). Rapid convergence is 
observed in the eigenenergies as we increase the size of the matrix in solving the equation 



(|10|). We choose, in the rest of this paper, 500 by 500 as a relevant size of the matrix in 
actual numerical calculations. In this case, the discrepancy in eigenenergies obtained by the 
shooting and finite difference methods is not greater than O.lmeV. 

3. RESULTS AND DISCUSSION 



The subband energies of an electron and a hole in GaAs/AlGaAs SQW are evaluated and 
the effects of carrier interaction and electric field are examined. The physical quantities used 
in obtaining numerical results are listed in table I. In this paper, we ignore the variation of 
the effective masses in the well and barrier regions of the SQW. And we consider the case 
that the bulk impurity density A^^ = 10^^/cm^ in the rest of the paper. 

3.1. Self Consistency Effects 

The potential profiles of the modulation-doped SQW in the presence of an external d. c. 
electric field F = 30kV/cm are shown in figure 3 for three different surface carrier densities 
Ng. In the figure, the vertical axis, especially the energy gap in the well is not drawn in 



scale. In the well region, the bottom of the conduction band moves upwards as one include 
the Hartree potential Vh{z), i.e., the interaction of a given electron with other electrons in 
the quantum well and also with the ionized impurities doped in the barrier region. Hence, 
the electron subband energies are also expected to increase due to the Hartree correction. 
In the valence band, however, the holes are more confined by the Hartree potential than 
electrons in the conduction band. This effect could cause a reduction of the hole subband 
energies. On the other hand, in modulation-doped SQW structures, there occurs a band- 
bending effect due to the dipoles formed between the plus (ionized donors) and the minus 
(electrons) charges and this effect would modify the subband energies. The effects of doped 
ionized impurities on the electron and hole subbands are different because of the opposite 
sign of the electric charges of an electron and a hole. 

Figure 4 shows electron subband energies E^n as a function of the carrier density Ng for 
four different values of the electric field strength F. The ground (n = 1) and first excited 
{n = 2) subband bottom energies are displayed, respectively, in the figure 4(a) and (b). The 
insets indicate that the electron subband energies are measured upwards from the bottom 
of the conduction band of the well with F = and A^^ = 0, which is taken the same as 
the midpoint of the conduction band bottom of the well in the presence of the field within 
a single-particle scheme; with Ns = . The Een decreases as one increases the magnitude 
of the electric field F. The reduction of the subband energies under the electric field is in 



agreement with observations of quantum confined Stark effect ||2^ . 

The Ng dependence of the subband energies is different in magnitude for different sub- 
bands. The subband energy E^i increases monotonically as A^^ increases, [figure 4(a)] 
However, the first excited subband energy Ee2 shows different behavior. [figure 4(b)] The Ee2 
increases initially as A^^ increases. For a given value of F, beyond a certain density A^*'^, -Ee2 
decreases until the subband n = 2 becomes quasi-bound to tunnel out the barrier located at 
the right-hand side. The ground subband energy Eei increases more rapidly than that of the 
first excited subband Ee2- This different behavior is related to the effect of band-bending 
on the barrier region. Figure 5 shows the probability density |'?/'e(-2)P of the ground(n = 1) 
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and first excited {n = 2) electron subbands for F = 30kV/cm and Ns = 0.9 x 10^^/cm^ and 
1.7 X 10^^/cm^. In tlie well region, as the surface carrier density A^^ increases, the magni- 
tudes of both |?/^e(^)P for both n = 1 and 2 decrease. On the other hand, the probability 
densities increase in the barrier region as the A^^ increases. The reduction of the |'?/'e(2)P in 
the well region results from the fact that the Hartree potential tends to locate the electrons 
near the interfaces, raising the corresponding confinement energies. However, the effect of 
band bending reduces the effective width of the potential barrier [see figure 3] and causes 
penetration of wave functions into the barrier region. The wave function penetration is more 
pronounced for the first excited subband than the ground subband. For the first excited 
subband, most of the electrons remain near the edges of the well for F < 50kV/cm even in 
the absence of the Hartree potential. But, the electrons in the ground subband, if F is not 
extremely large, such as, F < 150kV/cm, are well localized inside the quantum well by the 
confining potential Vt{z), in which the Hartree potential Vh{z) is excluded. 

Figure 6 shows the energy separation -Ee2-ei of the ground and first excited electron 
subbands as a function of A'^, for F = 0, 30, 50, and 70 kV/cm. The separation of electron 
subband energies decreases as A^^ increases because the effect of band bending is more 
pronounced in the first excited subband than in the ground subband. This result differs 
from the cases of Si inversion layers or GaAs single heterostructures p4|. In the cases of 



Si inversion layers and GaAs heterostructures, subband energy separations are observed to 
increase with increasing surface carrier density. This difference is conjectured to be related 
with the different boundary conditions and the effective potentials Vt{z). For the case of Si 
inversion layers, the potential barrier at the interface of Si and Si02 is well approximated to 
be infinitely high so that the wave functions vanish at the interface. For the case of single 
heteroj unction interface, the barrier height is finite. The potential functions Vt{z) in both 
cases are of triangular shape and infinite number of subbands are expected to be bound. 
However, in the SQWs, two heteroj unction boundaries exist and the shape of potential 
function is rectangular one with finite depth, hence, allowing only a few subbands to be 
bound. 
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The ground hole subband energies Ehi and first excited hole subband energies Eh2 cire 
shown, respectively, in figure 7(a) and (b) as a function of A^^ for various electric fields F. 
The hole subband energies are measured down from the midpoint of the valence band top 
of the well region with Ng = 0, as is indicated in the insets of figure 7(a) and (b). From the 
carrier-density dependence of the potential energy profile of the valence band [see figure 3] , 
the confinement of the holes is expected to increase as Ns increases. As one increases the 
carrier concentration of the SQW, the top of the valence band moves up so that the holes are 
more tightly bound. Hence, hole subband energies decrease with increasing surface carrier 
density. The electric-field dependences of Ehi and Eh2 are opposite to each other. The first 
excited hole subband is less bound for higher electric fields[see figure 7(b)]. 

Figure 8 shows the probability densities liphiz)]"^ of the ground(n = 1) and first 
excited(n = 2) hole subbands in the presence of an electric field F = 30kV/cm for two 
different surface carrier densities Ns = 0.9 x 10^^/cm^ and 1.7 x 10^^/cm^. We observe that 
the wave functions are more confined near the center of the well for higher A^^ and this 
confinement is more pronounced for the case of the ground subband. The first excited hole 
subband energies increase as F increases as is shown in figure 7(b). This peculiar behavior 
is caused by the fact that the first excited wave function is more confined to the center of the 
well as the electric field increases |^5|]. Figure 9 shows the energy separation of the ground 
and first excited hole subband energies Eh2-hi as a function of A^^ for four different values of 
the electric field strength F. We observe that the hole subband energies separation increases 
as one increases either Ng or F. This observation differs from the case of electron subbands 
[see figure 6] because of the different potential profiles[see figure 3]. In the zero-field case, 
the first excited hole subband energy (^ 33meV) is relatively small compare to the valence 
band offset (~ 160meV). Therefore, the effect of the Hartree potential on the first excited 
hole subband is much less than that on the corresponding electron subband. 
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3.2. Many Body Effects 



We include the effects of exchange and correlations in our calculation by employing a simple 
local exchange-correlation potential energy V^dz) given by equation (^. Strickly speaking, 
the density-functional formulation is good only for the calculations of the ground state 
energy and electron density distribution [rU-jT^. However, it has been used, in a good 
approximation, for the subband structure calculation The potential energy profiles of 
the conduction and valence bands of a SQW are shown, respectively, in figure 11(a) and (b) 
in the presence of an external d. c. electric field with (solid line) and without (dotted line) 
the exchange-correlation potential. The exchange-correlation potential slightly reduces the 
Coulomb repulsive force and electrons are pushed further toward the interface [|l^ . However, 
in this work, we observe that Vxc{z) lowers the conduction and valence band profiles of the 
SQW by negligiblely small amount. The magnitude of the Vxc{z) is so small, compared to 
other terms in Vt, that its effect on the total potential energy profiles is quite small. This 
result differs from the cases of Si inversion and accumulation layers [0,0] . 

The subband energy separations of Ee2-ei and Eh2-hi are shown, respectively, in figure 
11 (a) and (b) as a function of the surface carrier density Ng with (solid line) and without 
(dotted line) the exchange-correlation potential. The effect of the exchange-correlation 
potential is not significant and is only to increase subband energy separations slightly as is 
shown in the insets of figure 11(a) and (b). This increase comes from the fact the ground 
subband energy is lowered slightly more than that of the excited subbands, in agreement 
with the cases of Si layers The negligible effect of the exchange-correlation potential 
on the subband structure of the SQW system is conjectured to the relatively large (kinetic) 
confinement energy of the system, as compared with the cases of the single heterostructures 



1 2^ or Si space charge layers p3 |. 
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4. CONCLUSIONS 



We investigated the electronic properties of symmetrically-doped GaAs/AlGaAs single 
quantum wells in the presence of an external d. c. electric field perpendicular to the 
system and examined the effects of the field and carrier interactions on the system. 

We summarize the results as follows. First, we evaluated the subband energies and wave 
functions as a function of surface carrier density. A self-consistent calculation within a 
Hartree approximation shows that, as the surface carrier density increases, the ground elec- 
tron subband energy increases monotonically, but the first excited subband energy starts to 
decrease beyond a certain density. However, the energy separation of the lowest two electron 
subbands is reduced as the surface carrier density increases. For the hole subbands, as the 
surface carrier density increases, individual subband energy decreases, but their energy sep- 
aration increases. This different behavior in the carrier-density dependence of the electron 
and hole subbands results from the different confinement effects for two different charge car- 
riers. Second, the exchange-correlation effect is examined within a local density-functional 
approximation. A small reduction of the subband energies and an slight enhancement of 
the subband energy separation is observed. However, the many-body effect is found not so 
significant. This result is different from the cases of the Si inversion layers or GaAs single 
heter ostructures . 
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Figure Captions 

Figure 1. Schematic diagram of an idealized quantum well structure of width L and barrier 
height Vo in the presence of an external d. c. electric field F. The wave functions are to 
vanish at —Lb and +Lb for the convenience of a calculation. The width of symmetrically 
doped region is indicated by db in the barrier regions. 

Figure 2. The variation of the electron subband energies as a function of the size of matrix. 
The results are the case of undoped SQW in the absence of the electric field. The width of 
the well is 85A. 

Figure 3. The potential energy profiles of modulation-doped SQW in the presence of an 
electric field F = 30kV/cm for different surface carrier densities Ng. The width of the well 
is 85A. The vertical axis, especially, the energy gap in the well region is not drawn in scale. 

Figure 4. The electron subband energies as a function of the surface carrier density A^^ 
for various electric fields F. The inset indicates that the energy is measured up from the 
midpoint of the conduction band bottom in a single-particle scheme in the presence of the 
fields, (a), the ground electron subband. (b). the first excited electron subband. 

Figure 5. The probabihty density of the ground and first excited electron subbands under 
an electric field F — 30kV/cm for two different surface carrier concentrations Ng. The width 
of the well is 85 A. 

Figure 6. The energy separation of the ground and first excited electron subbands as a 
function of the surface carrier density A'^g for four different values of electric field strength. 
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Figure 7. The hole subband energies as a function of the surface carrier density Ns for var- 
ious electric fields F. The inset indicates that the hole subband energies are measured down 
from the midpoint of the valence band top in the well region in a single-particle scheme in 
the presence of the fields, (a), the ground hole subband. (b). the first excited hole subband. 

Figure 8. The probability density of the ground and first excited hole subbands under an 
electric field F — 30kV/cm for two different surface carrier concentrations Ng. The width 
of the well is 85 A. 

Figure 9. The energy separation of the ground and first excited hole subbands as a function 
of the surface carrier density Ng for four different values of electric field strength F. 

Figure 10. The potential energy profiles with (soild fine) and without (dotted fine) the 
exchange-correlation potential in the presence of an electric field F — 30kV/cm for a con- 
stant surface carrier density Ng = 1.7 x 10^^/cm^. The width of the well is 85A. (a), the 
bottom of the conduction band. (b). the top of the valence band. 

Figure 11. The energy separation of the lowest two subbands with (solid line) and without 
(dotted line) the exchange-correlation potential as a function of the surface carrier density 
Ng for an electric field F — 30kV/cm. The solid and dotted lines are superposed nearly 
on top of each other, (a), the electron subband energy separation, (b). the hole subband 
energy separation. 
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TABLES 

TABLE I. Physical quantities used in numerical calculation 



(mo : free electron mass) 



Quantum Well Material 


GaAs/AlGaAs SQW 


Well width 


85A 


Electron effective mass 


0.067mo 


Hole effective mass 


0.45mo 


Conduction band offset 


240meV 


Valence band offset 


leOmeV 


Background dielectric constant 


13.1 


Doped impurity density, Nd 


1.0 X 10^^ /cm^ 
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